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Abstract 

A generic synchrotron external shock model is the widely preferred 
paradigm used to interpret the broad-band afterglow data of gamma- 
ray bursts (GRBs), including predicted observable signatures from a 
reverse shock which have been confirmed by observations. Investiga¬ 
tions of the nature of the reverse shock emission can provide valuable 
insights into the intrinsic properties of the GRB ejecta. Here we briefly 
review the standard and the extended models of the reverse shock emis¬ 
sion, discussing the connection between the theory and observations, 
including the implications of the latest observational advances. 


1 Introduction 

Gamma-ray bursts (GRBs), which are the most extreme explosive events 
in the universe, generally present two phenomenological emission phases: 
an initial prompt y-ray emission and a longer-lived broadband afterglow 
emission. Regardless of the nature of the progenitor and the central engine, 
the radiation of the GRBs is believed to be caused by the dissipation of 
the kinetic energy of a relativistic jet which is beamed towards Earth (for 
reviews, see Ref EMEUS]). Although the detailed physics of the prompt 
y-ray emission is still uncertain, mainly owing to the poorly understanding 
composition of the GRB jet (e.g., the degree of magnetization) [6], a generic 
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synchrotron external shock model is the most widely accepted paradigm for 
interpreting the broad-band afterglow data M Ena eh E2]. 

The external shock model is based on a relativistic blastwave theory that 
describes the interaction between the GRB jet (i.e. the ejecta) and the cir- 
cumburst medium (for detailed reviews, see |T3]). During the interaction, 
two shocks naturally develop. A long-lived forward shock sweeps up the am¬ 
bient medium, which gives rise to the long-term broad band afterglow; and 
a short-lived reverse shock propagates into the GRB ejecta, which can give 
rise to a short-term optical/IR flash and a radio flare. In the pie-Swift era, 
the forward shock signal was found to successfully represent a large array of 
late-time afterglow data mmmmmmmm although moderate 
revisions are sometimes required m for the more complicated afterglow be¬ 
haviors [22 | l23 | IM1 l25l [26] . After the launch of NASA’s dedicated GRB mis¬ 
sion Swift m, unprecedented new information about GRB afterglows was 
revealed [23 E9] EI2 EH 1MI ES] , especially in the early phases, thanks mainly 
to the rapid slewing and precise localization capability of its on-board X-Ray 
Telescope (XRT) [M). It was found that a number of physical processes are 
needed to shape the observed lightcurves [Ml HU], including, e.g., the sugges¬ 
tion that the X-ray afterglow is a superposition of the conventional external 
shock component and a radiation component that is related to the late cen¬ 
tral engine activity [H [29l [301 ESI (Ml EH [38l (Ml SOI HD HU [Ml HU SSI S6] . 
In any case, the external forward shock still remains the basic theoretical 
framework to interpret the broad band afterglow signals. It is elegant in 
its simplicity, since it invokes a limit number of model parameters (e.g. the 
total energy of the system, the ambient density and its profile), and has 
well defined predicted spectral and temporal properties. However, it lacks 
the capability to study some detailed features of the GRB ejecta, such as 
the composition, since its radiation comes from the shocked medium rather 
than the ejecta materials. 

The reverse shock, on the other hand, should heat the GRB ejecta within 
a short period of time, contributing another important aspect to the external 
shock emission signature. The hydrodynamics of reverse shock propagation 
in a matter-dominated shell and its corresponding radiation features were 
studied in great detail [91 QHl EH m\ prior to the expected signals being 
discovered. In the pr e-Swift era, some cases with very early optical flashes 
(e.g. GRB 990123 [22]; GRB 021004 [25]; GRB 021211 [Ml [19]) or early 
radio flares [50] were detected, which generally agreed well with the predicted 
reverse shock emission [511E21ESI [Ml ESI ESI EH [SHI EHl HEH EH E21 [63]- 

However, there are also some observations which challenge the simple 
reverse shock prediction. For instance, the early optical emission of GRB 
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030418 [BI] does not agree with the predicted reverse shock behavior; fur¬ 
thermore, rapid optical follow-up observations for some bursts reveal the 
so-called “optical flash problem”, e.g. upper limits of 15 mag were estab¬ 
lished for specific observed bursts, instead of detecting the expected reverse 
shock emission [HS, §3 E3 E5]- In order to better interpret the observa¬ 
tional results, the simple reverse shock model was extended to accommodate 
more realistic conditions than what was initially assumed. E.g., the ambient 
medium might be a stellar wind (or in general have a profile n oc r~ k ) rather 
than being a uniform interstellar medium [69, i 70l fTTl 172] ; the reverse shock 
propagation speed might be semi-relativistic instead of ultra-relativistic or 
non-relativistic [62]; the GRB ejecta might be magnetized, which could en¬ 
hance the signal when the magnetization is moderate, or completely suppress 
the signal when magnetization degree is large enough [731E3 ES ESI El [78]; 
the GRB outflow may carry a good fraction of electron-position pairs or 
neutrons which could alter the early afterglow behavior m [80]; consid¬ 
ering a more complicated stratification profile of the ejecta, e.g., with a 
nonuniform Lorentz factor, luminosity and density, the reverse shock emis¬ 
sion could have a richer set of features, including being able to reproduce 
the canonical X-ray lightcurves as observed by Swift as long as the forward 
shock emission is suppressed EDIE3E3EDE3 [86] . Besides these model 
modifications, some new signatures for reverse shock were also proposed, 
such as sub-GeV photon flashes and high-energy neutrino emission [87l [88 ], 
early X-ray and gamma-ray emission from synchrotron self-Compton (SSC) 
in the reverse shock region or cross inverse Compton (IC) between the elec¬ 
trons and photons from the forward shock and reverse shock |89[ 190] , or a 
polarization signature that offers the possibility to diagnose the structure of 
the magnetic fields in the GRB ejecta, etc. 

Before the launch of Swift , the observational data was not ample or 
detailed enough to comprehensively test these reverse shock models nor to 
study the ejecta properties through the reverse shock signatures. A good 
sample of early afterglow lightcurves which would allow a detailed study 
of GRB reverse shocks was one of the expectations from the Swift mission 
mm- After ten years of successful operation of Swift , it is now of great 
interest to revisit this problem and to see how much progress has been made. 

The structure of this review is as follows: we first summarize the models 
for the reverse shock emission, including the standard synchrotron external 
shock model in section 2, and discuss the extended models in section 3. In 
section 4 we illustrate how to identify in practice the reverse shock signals 
present in the observational data, and how to use such signals to study 
the GRB ejecta properties. The current observational results and their 
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implications are collected in section 5. We conclude with a brief discussion 
of the prospects for future reverse shock studies. 

2 Standard modeling of the reverse shock emission 

2.1 Model description 

Consider a uniform relativistic coasting shell with rest mass Mo, energy 
E, initial Lorentz factor rj = E/Mqc 2 , and observed width A, expanding 
into the circumburst medium (CBM) described by a density profile n{r) = 
Ar~ k , 0 < k < 4. A pair of shocks will develop, namely, a forward shock 
propagating into the medium and a reverse shock propagating into the shell. 
The two shocks and the contact discontinuity separate the system into four 
regions: (1) the unshocked CBM (called region 1 hereafter), (2) the shocked 
CBM (region 2), (3) the shocked shell (region 3), and (4) the unshocked 
shell (region 4). Synchrotron emission is expected from regions 2 and 3, 
since electrons are accelerated at the shock fronts via the lst-order Fermi 
acceleration mechanism and magnetic fields are believed to be generated 
behind the shocks due to plasma instabilities (for forward shock) [91] or 
shock compression amplification of the magnetic field carried by the central 
engine (for reverse shock). 

An evaluation of the hydrodynamical and thermodynamical quantities 
for the region 2 and 3, namely, 'ji,rii,pi and e* (bulk Lorentz factor, particle 
number density, pressure and internal energy density, with i denoting the 
region number), allows one to straightforwardly calculate the instantaneous 
synchrotron spectrum at a given epoch, as well as the flux evolution in 
time (the lightcurve) for a given observed frequency. In doing this, it is 
customary to introduce parametrizations for the microscopic processes, such 
as the fractions of the shock energy that go into the electrons and into 
magnetic fields (e e and e#), and the electron spectral index ( p ). Ref [13] 
gives detailed examples about such calculations and provides a complete 
reference for all the analytical synchrotron external shock afterglow models 
by deriving the temporal and spectral indices of all the models in all spectral 
regimes. In order to review the reverse shock related features, we give here 
a brief summary of the dynamical properties of region 3 for various models. 

In general, region 3 will evolve through two different phases, i.e., before 
the reverse shock crossing the shell (at T x ), and after the reverse shock 
crossing. The dynamical solution depends on the relativistic nature of the 
reverse shock, which can be characterized by the dimensionless parameter 
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£ = (1/ A) 1//2 r? ( 4 03 72], where l = ( 4 nAm % ) 3 k ^ ie Sedov 

length (at which the swept-up medium’s rest-mass energy equals the initial 
energy E of the shell). If £ <C 1, the reverse shock is ultra-relativistic 
(thick shell regime), while if £ 1, the reverse shock is Newtonian (thin 
shell regime). Between these two extreme limits, the reverse shock can be 
considered semi-relativistic when £ is of the order of unity [62J. Combined 
with the different (generic) types of CBM, i.e., constant density interstellar 
medium (ISM) model [k = 0), stellar wind model ( k = 2) and general 
stratified wind model (0 < k < 4), seven different regimes have been studied 
in the literature [47l 153] l62l [Hill 1701 [7T1 [72] . These are: 1) thick shell ISM 
(£ <C 1, n\ oc r°); 2) thin shell ISM (£ S> 1, n\ oc r°); 3) thick shell stellar 
wind (£ < 1, nj oc r~ 2 ); 4) thin shell stellar wind (£ 1, n\ oc r~ 2 ); 5) 

thick shell general stratified wind (£ <C 1, n\ oc r~ k ); 6) thin shell general 
stratified wind (£ 1, n\ oc r~ k )\ 7) semi-relativistic reverse shock ISM 

(£ ~ 1, n\ oc ?’°). Below we summarize the results in the literature for these 
different regimes. 

1) Thick shell ISM (£ <C 1, m oc r°) 07JO3J 

In this case, the reverse shock crossing time can be estimated as T x = 
A/c, which is independent of the CBM (applied to all thick shell regimes 
below). Before T x , the dynamic variables of region 3 in terms of the observer 
time t = r/ 2 c 7 | are 
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where Nq = E/r]m p c 2 is the total number of electrons in the shell. Since 
the shocked regions (region 2 and 3) should be extremely hot, the energy 
density term is degenerate with the pressure term as e% = 3p%. 

After T x . the profile of the shocked medium in region 2 begins to ap¬ 
proach the Blandford-McKee (BM) self-similar solution (921193]. Since re¬ 
gion 3 is located not too far behind region 2, it should roughly fit the BM 
solution, which is verified numerically as long as the relativistic reverse shock 
can heat the shell to a relativistic temperature [94] . The BM scaling thus 
can be applied to the evolution of the shocked shell, 
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/ t \ -13/12 

P3 = P3(T x ) l — j , N 3 = No. 


( 2 ) 


Note that the number of the shocked electrons is constant after the shock 
crossing since no electrons are newly shocked. 

2) Thin shell ISM (£ /§> 1, n\ oc r°) [47, 5l] 

In a thin shell case, the reverse shock is too weak to decelerate the 
shell effectively. T x can be estimated by the deceleration time of the ejecta 
(applied to all thin shell regimes below) 


7-e — tdec 
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Before T x , the scalings for the dynamic variables of region 3 is given by 
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After T x , the Lorentz factor of the shocked shell may be assumed to 
have a general power-law decay behavior 73 oc r~ 9 [52l ;53j. The dynamical 
behavior in region 3 may be expressed through the scaling-laws 


73 OC t -3/(l+2 9 ) )n30ct -6(3 +9 )/7(l + 2 p)) 

p 3 OC i -8( 3 +s)/n i +2 P ) )7 . oc t 1 /(i+ 2 fl) j JV e)3 oc t°. (5) 


For the ISM case, numerical studes showed that the scalings with g ~ 2 fits 
the evolution [53], e.g., 

73 = 73( T -)(^) 2/5 , n 3 = n 3 (T*)Qr) ^, 
f t \ " 8/7 

P3 = P3(T x ) , N 3 = N 0 . ( 6 ) 


3) Thick shell stellar wind (£ <C 1, n\ oc r 2 ) [69, !7D] 
Similar to regime 1, before T x , we have 
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After T x , assuming a BM self-similar adiabatic solution for the evolution 
of the shocked shell m , the relevant hydrodynamic variables are given by 

73 = 73 (Tx) (^) ' ,n 3 = n 3 (T x )(j^j ■ 

( t \ ~ 3/2 

P3 = P3{T x ) J , iV 3 = iVo- (8) 

4) Thin shell stellar wind (£ S> 1, n\ oc r -2 ) [71] 

In this case, the evolution of the hydrodynamic variables before the time 
T x are 


73 = V, ri3 = 


7Ar/ 6 ( t 
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After the reverse shock crosses the shell, the scaling law for regime 2 still 
applies, except g = 1, namely 

73 = 73 (T x )(pj 1/3 , n 3 = n 3 (T x )(^ ^ 


P3 = P3(T x ) 


t \ -32/21 


, N = N 0 . 


( 10 ) 


5) Thick shell general stratified wind (£ <C 1, ni oc r~ k ) 

Before the reverse shock crosses the shell, the hydrodynamical evolution 
of the reverse shock can be characterized by 


t \ —(2—*)/2(4—fc) / t \ -(6+fc)/2(4 -k) 
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Pz{T x ) = 3^{T x )n 3 (T x )m p c 2 . 


After the reverse shock crosses the shell, again with a BM self-similar 

2fc —7 4fc —26 2fc —13 o 

solution, one gets 73 oc r 2 , p 3 oc r 3 , n 3 oc r 2 , and t oc r/ygc. Thus, 

the hydrodynamic evolution of the reverse shock after crossing the shell is 
characterized by 


(12) 


/ t \ (2fc-7)/4(4-fc) / t \ (2fc—13)/4(4—fc) 

73 = 13 {T X ) 1 — J , n 3 = n 3 (T x ) f — J 

/ * \ (2fc—13)/3(4—fc) 

p 3 = p 3 (r x ) f — J , iv 3 = iv 0 . (13) 

6 ) Thin shell general stratified wind (£ 1, ni oc r~ k ) [72] 

In this case, before T x , the hydrodynamic evolution of the reverse shock 
can be characterized by 


where 
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= 3 (734 (T x ) - 1) n 3 (T x )m p c 2 


(14) 


(15) 


After the reverse shock crosses the shell, the scaling law for regime 2 
should be still relevant, except that the value of g has not been studied in 
detail. 

7) Mild relativistic reverse shock ISM (£ ~ 1, m oc r°) [62] 

In this case, a simple analytical solution is no longer achievable. The 
nature of the reverse shock is determined by £ and another parameter a, 
which is the ratio of the Lorentz factor of the shocked matter to r/, 


« = 73 / 7 - 


(16) 
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Here a can be derived directly from the relativistic jump conditions m ■ 


(12/£ 3 - l)a 4 + 0.5a 3 + a 2 + 0.5a -1 = 0. 


(17) 


The reverse shock reaches the back of the shell at 


T x = — (1 + 0.5A4£ 3/2 ) 

c 


(18) 


where Aft = 1.4 is a numerical correction factor to the analytic estimates 
[62 j. At this stage, 



When t < T x , the dynamical variables of region 3 can be determined by 
parameterizing all the quantities according to the fraction of the reverse 
shock crossing the shell, /: 


A(/) oc E(f) oc /, C(/) oc / _1/3 , 
Hf) oc f l/2 Af) oc /(I + 0.5A/)£(/) 3/2 ). 


( 20 ) 


At t > T x , the hydrodynamical evolution becomes almost independent of 
e m , therefore the solutions for the dynamic variables of region 3 become 
the same as in regime 2. 

2.2 Emission evolution 

The instantaneous synchrotron spectrum at a given epoch can be described 
with three characteristic frequencies v a (self-absorption frequency), v m , and 
v c (the cooling frequency), and the peak synchrotron flux density Tj, iinax [11 j • 
Based on the dynamical solution for specific situations, one can calculate the 
temporal evolution of these characteristic parameters and then determine 
the flux evolution in time (the lightcurve) for a given observed frequency. 
Since the reverse shock emission is expected to be prominent in the optical 
band at early stage, here we give a brief description for the morphology of 
early optical afterglow lightcurves. 

It is shown that for reasonable parameter spaces, shortly after (or even 
during) the prompt emission phase, both forward shock and reverse shock 
emission would enter into the “slow cooling” regime {v c < v m ) [HI Ell¬ 
in the following, we will take slow cooling for both reverse and forward 
shock emission, so that the shape of the lightcurve essentially depends on 
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the relation between Vm and n op t, where the superscript r and / represent 
reverse and forward shock respectively. 

For thin shell case, the evolution of I'm reads 

z 'L oc t° (t < T x ), oc t~ 3/2 (t > T x ), 

U r m oc t 6 (t < T x ), z4 oc t- 54 / 35 (t > T x ). (21) 

As shown in Figured^,, when i'll (T x ) is larger than u opt , I'm would cross 
the optical band once for the forward shock (at tf) and twice for the reverse 
shock (at t T} i and t r p). In this case, we have (shown in Figure [T)o) 

Ft <x t 3 (t < T x ), Ft oc t 1/2 (T x < t < t f ), Ft (x t" 3 ^" 1 )/ 4 (t > t f ), (22) 

and 

FI oc t( 6 P“ 3 )/ 2 (t < t r , i), FI oc t~ 1/2 (t r> i < t < tr }2 ), FI oc ^( 27 p+ 7 )/ 35 (t > i r>2 ).(23) 

When Vni {T x ) is smaller than z/ opt , there is no v m crossing and the lightcurves 
for both shocks peak at T x . In this case, we have 

Fl (x t 3 (t < T x ), Ft (x t- 3 ^" 1 )/ 4 (t > T x ), 

FI oc t (6p ~ 3)/2 ( t < T x ), Fl oc H 27 p+ 7 )/ 35 {t > T x ). (24) 

Depending on their shapes and relative relations between the forward shock 
and reverse shock emission, the early optical light curves could be distributed 
into different morphological types, we will discuss this in detail in section 

HU 

For thick shell case, the evolution of ul/ reads (shown in Figured}:) 

u f m oc t _1 (t < T x ), v f m oc t~ 3/2 (t > T x ), 

< oc t° (t < T x ), v r m (x f 73 / 48 (t>T x ). (25) 

When I'm (T x ) is larger than n opt , I'm would cross the optical band once for 
both forward shock (at t/p) and reverse shock (at t r ). In this case, we have 
(shown in Figured}!) 

Fl oc f 4/3 (t < T x ), Fl oc t 1/2 (T x <t< tf), Fl oc t~ 3{ - p ~ l) F (t>t f ), (26) 
and 

Fl oc t 1/2 (t < T x ), Fl oc t" 17/36 (T x < t < t r ), Fl oc f( 73 P+ 21 )/96 > t r ).(27) 
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When I'm (T x ) is smaller than I'optj there is no u m crossing for reverse shock 
but one time crossing for forward shock (at tf 2 )- In this case, we have 

H (x t 4 / 3 (t < t ft2 ), Ft oc (i />2 < t < T x ), Fl oc f 3 ^- 1 )/ 4 (t > T x ),( 28) 

and 


F r v oc t 1/2 (t < T x ), F r v oc t-( 73 P+ 21 )/ 96 (* > T x ). (29) 



Thick shell case 



Figure 1: Illustration of the v m evolution (left panels) and optical lightcurves 
(right panels) for both forward shock (blue lines) and reverse shock (black 
lines) emission, from Ref [95]. Red circles on lightcurve indicate the points 
for comparison in order to categorize the lightcurve types [95]. 
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3 Extended models of the reverse shock emission 


3.1 Reverse shock emission from magnetized GRB ejecta 

It has been suggested that the GRB ejecta is likely to be magnetized (see 
Ref [5j for a recent review). Although the degree of magnetization is still 
unknown, it is usually quantified by through the parameter cr, the ratio of 
the electromagnetic energy flux to the kinetic energy flux. The existence of 
magnetic fields in the ejecta will influence at least two aspects of the reverse 
shock characteristics, i.e., the hydrodynamical solutions for the shocked shell 
region and the reverse shock emission level. 

Under ideal MHD conditions and with a more accurate approach to ac¬ 
count for the modifications in the shock jump conditions when magnetic 
fields are involved, a rigorous analytical solution for the relativistic 90° 
shocks was carried out and several interesting conclusions were suggested 

m- 

• A strong reverse shock still exists in the high-cr regime, as long as the 
shock is relativistic. For typical GRB parameters, the reverse shock 
could form when cr is as high as several tens or even hundreds, which is 
supported numerically by solving the one-dimensional Riemann prob¬ 
lem for the deceleration of an arbitrarily magnetized relativistic flow 

• The dynamical evolution of region 3 can be still categorized into the 
thick and thin shell regimes, except that the pivotal parameter to 
separate the two regimes now becomes a. At larger cr-value, the thick 
shell regime greatly shrinks and the reverse shock emission peak is 
broadened in the thin shell regime due to the separation of the shock 
crossing radius and the deceleration radius. Such novel features could 
be useful for diagnosing the magnetization degree of GRB ejecta. 

• The reverse shock emission level should initially increase rapidly as cr 
increases from below, until reaching a peak around cr ~ 0.1 — 1, and 
decreases steadily when a > 1. The decrease of the emission level is 
caused not only because the reverse shock becomes weaker, but also 
because the total kinetic energy fraction in the flow gets smaller. Sep¬ 
arate investigations of the reverse shock emission powered by mildly 
magnetized (cr ~ 0.05 — 1) GRB ejecta were also carried out numeri¬ 
cally [96] . and similar results were achieved. In that work [96], both 
ISM and stellar wind CBMs were considered, and it turns out that be¬ 
fore the reverse shock crosses the ejecta, the relevant R-band emission 
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flux increases rapidly for the ISM medium case, but for the wind case 
it increases only slightly, which is similar to non-magnetized scenario. 
Recently, multi-band GRB afterglow lightcurves for magnetized ejecta 
have been calculated with high-resolution relativistic MHD simulations 
coupled with a radiative transfer code cans], and it is suggested that 
for typical parameters of the ejecta, the emission from the reverse 
shock peaks at magnetization values a ~ 0.01 — 0.1 of the flow, and 
that it is greatly suppressed for higher c-values. 

• In the high o'-value regime, a sufficient magnetic energy has not yet 
been transferred to the ISM at the end of the reverse shock crossing, 
since the magnetic pressure behind the contact discontinuity balances 
the thermal pressure in the forward shock crossing. The leftover mag¬ 
netic energy would eventually be injected into the blastwave or dissi¬ 
pate into radiation at some point and provide additional signatures to 
the afterglow lightcurve [731176]. 

3.2 Reverse shock emission from pair-rich or neutron-fed 
GRB ejecta 

Besides magnetic fields, other components of the GRB ejecta, if present, 
could also alter the reverse shock emission features, such as e± pairs and 
neutrons [26] [5 8 ]. 

The intrinsic GRB spectrum may extend to very high energy, so that 
the optical depth to 7 — 7 absorption for the most energetic photons at the 
high-energy end of the spectrum may exceed unity. In this case, intense pair 
production may occur in the prompt emission phase and pairs remain 
in the fireball, with the same bulk Lorentz factor as the fireball (static 
in the comoving frame). Since the pair will also share energy in the 
reverse shock, the reverse shock emission spectrum is altered, and the peak is 
softened to lower frequencies. It turns out that a pair-rich reverse shock gives 
rise to stronger radiation in the IR band, instead of the optical/UV emission 
in the case where pair-loading is negligible [26] . The optical afterglow signal 
may suffer significant dust obscuration since long GRBs are usually expected 
to occur within star forming regions; observable IR flashes could test this 
issue, provided IR detector can be slewed rapidly enough to respond the 
GRB trigger [2Bj. 

It has also been pointed out that GRB ejecta may contain a significant 
fraction of neutrons pH] EH EH 199] . which would cause much more com¬ 
plex dynamics for the system than in the neutron-free case. In general, 
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the neutron shells (IV-ejecta) would freely penetrate through the charged 
ion shells (/-ejecta) in front of them, and would separate from the /-ejecta 
more and more, while the /-ejecta suffer deceleration from internal shocks. 
The IV-ejecta would decelerate by collecting ambient medium and the mass 
of fast neutrons would decrease as the result of /3-decay. The neutron decay 
products and the shocked medium will form a new ejecta (T-ejecta) that 
follows behind the IV-ejecta and the interactions between these three ejecta 
would give rise to rich radiation features. For an ISM type medium, the 
T-ejecta moves faster than the /-ejecta, so that the T-ejecta would first 
interact with the iV-ejecta or ambient medium, but the reverse shock emis¬ 
sion in this stage would be out-shined by the forward shock emission. Later 
on, the /-ejecta would catch up the T-ejecta and a prominent bump signa¬ 
ture around tens to hundreds of seconds would show up, which is mainly 
dominated by the refreshed reverse shock emission. For a stellar wind type 
medium, /-ejecta would pick up the T-ejecta first and then collide with the 
IV-ejecta and ambient medium. In this case, three components contribute 
to the final emission, i.e. the forward shock emission, the reverse shock 
emission from the shocked /-ejecta and the shocked T-ejecta emission. A 
typical neutron-rich wind-interaction lightcurve is characterized by a promi¬ 
nent early plateau lasting for ~ 100 s, followed by a normal power-law decay 

m- 

3.3 High energy photons and neutrinos from reverse shock 

Since the number of heated electrons in region 3 is 77 (10 2 — 10 3 ) times 
higher than in region 2, a strong synchrotron self-Compton (SSC) emission 
in region 3 is expected, especially when reverse shock emission is prominent 
Msnuino]. The SSC emission feature is essentially determined by the 
random Lorentz factors of the electrons r j e , since the seed photons mainly 
are concentrated in the optical band. When y e is of the order of 1000 or even 
higher, the SSC emission from the reverse shock could dominate over the 
synchrotron and other IC emissions in the energy bands from tens of MeV to 
tens of GeV, while the cross-IC (and/or the forward shock SSC emissions), 
becomes increasingly dominant at TeV energy bands [891 [90]. When y e is 
of order 100, if the SSC process dominates the cooling of shocked electrons, 
the majority of the shock energy would be radiated in the second-order 
scattering at 10 — 100 MeV, and the first-order scattering may give rise to 
X-ray flares in the very early afterglow phase [ lOO j. In this case, the optical 
flash (due to synchrotron) is highly suppressed. 

On the other hand, it has been proposed that when GRBs erupt in 
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a stellar wind, usually the region 2 and region 3 still have overlap with 
the prompt MeV 7 -ray emission site at the reverse shock crossing phase 
[80, 101]. Such overlapping could lead to significant modifications of the 
early afterglow emission, since the dominant cooling process for the electrons 
is likely to be the IC process with the MeV photons jlOl] . Due to the close 
overlap of the MeV photon flow and the shocked regions, the newly up- 
scattered high energy photons would be absorbed by the MeV photons to 
generate pairs, and then re-scatter the soft X-rays to power a detectable 
sub-GeV signal [80]. Other than that, 10 14 eV neutrino emission is also 
expected from interactions between shocked protons and the MeV photon 
flow [80]. Alternatively, high energy neutrinos are also expected from reverse 
shocks as the GRB jets crossing the stellar envelop, either for choked or 
successful relativistic jets [ 102 ] , 

3.4 Long lasting reverse shocks 

In the standard model, a uniform distribution of the bulk Lorentz factors 
in the GRB ejecta is assumed. However, in principle GRB ejecta could 
have a range of bulk Lorentz factors, so that the inner (lower 7 ) parts may 
carry most of the mass, or even most of the energy, e.g. 7 Me 2 oc 7 _s+1 
[83 , 103 , 1104] . In this case, the low Lorentz factor part of the ejecta will 
catch up with the high Lorentz factor part when the latter is decelerated 
by the ambient medium, thus a long-lasting weak reverse shock could de¬ 
velop, until the whole ejecta has been shocked. Analogously to the standard 
model, this process could also be classified analytically into two cases: the 
thick shell case and the thin shell case [83] , and it turns out in the thick shell 
case, the reverse shock is strong and may give rise to the plateau observed 
in the early optical and X-ray afterglows [83]. Considering more compli¬ 
cated stratification profiles for the ejecta properties (e.g., Lorentz factor, 
luminosity and density), the long lasting reverse shock emission could be 
endowed with a richer set of features, including reproducing the canonical 
X-ray lightcurve as observed by Swift , as long as the forward shock emission 
can somehow be suppressed [83] [85]. 

3.5 Polarization of reverse shock emission 

If the GRB ejecta contains large scale ordered magnetic fields, the prompt 
7 -ray emission and the reverse shock emission should be polarized m- 
However, aside from any instrumental difficulties, making unequivocal po¬ 
larization determinations that prove this are still challenging [ 1051 i 106] . Fur- 
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thermore, a high degree of linear polarization in the prompt 7 -rays is also 
possible in the presence of a random magnetic field, arguably originating 
in electromagnetic instabilities that develop at the collisionless shock [ 10*7 ]. 
In any case, polarization measurements of the reverse shock emission could 
place strong constraints on the strength, and perhaps also the structure of 
the magnetic field within the GRB outflow. The RINGO detector on the 
Liverpool Telescope has reported an optical polarization of GRB 090102 
(P = 101%) [ 108] and GRB 060418 (P < 8 %) [109], but a larger sample is 
definitely needed to give general discussion on the properties of GRB outflow 

M- 

4 Connection between theory and observations of 
the reverse shock emission 

4.1 Theory predictions of observational features of reverse 
shocks 

According to the standard external shock theory, reverse shocks would mainly 
contribute to the early optical afterglow (if not suppressed) [13]. For the 
ISM model, the early optical lightcurve of the reverse shock would increase 
proportional to t 5 (thin shell case) or t 1 / 2 (thick shell case), and then de¬ 
crease with a general slope ~ t ~ 2 [SUES]- For the wind model, the lightcurve 
would increase initially with slope f 5 / 2 when synchrotron self-absorption be¬ 
comes important in this case, and then rising with slope 1/2 for both thin 
and thick cases, to finally decrease with a slope ~ t” 3 , determined by the 
angular time delay effect m- 

The morphology of early optical afterglows essentially depend on the 
relative relation between the forward shock and reverse shock emission. In 
general, the early optical afterglows for constant density medium model were 
usually classified into three types (see Figure 2): 

• Type I: re-brightening. Two peaks emerge in this type of lightcurve. 
The first peak is dominated by the reverse shock emission, and the 
re-brightening signature comes from the forward shock emission. The 
temporal index for the re-brightening depends on the specific forward 
shock model and the spectral regimes, which are collected in Ref Em 

• Type II: flattening. In this case, the forward shock emission peak is 
under the reverse shock component, and the decaying part of the for¬ 
ward shock emission shows up later when the reverse shock component 
is getting fainter more rapidly. 
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• Type III: no reverse shock component. Two reasons may be respon¬ 
sible for this, one being that the reverse shock component is weak 
compared with the forward shock emission; the other being that the 
reverse shock component is completely suppressed for some reason as 
proposed by some extended models (see section ED, such as magnetic 
fields dominating the ejecta m, e ± pair effects [26j, or SSC process 
in the reverse shock region |521 EOj, lOOj . 

Recently, it is suggested that an insight on the r'm(T'z) value could lead to 
strong constraints on relevant afterglow parameters m, so that the forward 
shock dominated cases (Type III) should be redefined into two categories: 

• Type III: forward shock dominated lightcurves without v m crossing. 

• Type IV: forward shock dominated lightcurves with v m crossing, 

4.2 Identification of reverse shock emission from observa¬ 
tional data 

Based on the theoretically predicted features, once prompt optical observa¬ 
tions are obtained, the reverse shock components could be identified with 
the following procedure: 

1) Compare the first optical observation time t s and the 7 -ray duration 
Tqq. If t s < T 90 , check the variability level of the optical signal. For cases 
with significant variability, ascertain the relation between optical variability 
and 7 -ray variability with correlation cross checking method. Bursts with 
t s > T90, or t s < T90 but with weak variability (or with significant variability 
but no correlation with 7 -ray signal) may be taken as candidates for having 
a reverse shock signal. It is worth pointing out that variability within a 
certain level may be explained within the external shock framework, such 
as invoking density fluctuation, inhomogeneous jets, or neutron decay sig¬ 
natures, etc | 88 l llllj 11121ITI31 . Information from other observational bands 
(radio, X-ray, and high energy 7 -rays) would be helpful to make a stricter 
selection between cases. 

2) Fitting the optical lightcurve with a multi-segment broken power law 
function. If the initial decay slope of the signal is close to t ~ 2 (ISM) or 
t~ 3 (wind), check whether the following decay or rising slopes are consistent 
with the forward shock predictions m and classify the candidate bursts as 
one of the four types defined above. 

3) Plot the multi-band spectrum of the early afterglow, if possible, and 
verify if there is evidence for the existence of two components, e.g. forward 
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Figure 2: Theoretically expected early optical afterglow lightcurves from re¬ 
verse plus forward shock emission, and illustrative diagram of three classified 
types, (a) from Ref [58j; (b) from Ref [70]; (c) from Ref [110] . 


shock component (usually peaks at X-ray) and reverse shock component 
(usually peaks at optical). 

4.3 Constraints on theoretical parameters from observational 
results 

Valuable results may be expected in the case of bursts where multi-band 
(instead of only X-ray) early afterglow observations are available, especially 
for the properties of the GRB outflow itself. For cases with identifiable re¬ 
verse shock component, several important pieces of information, if available, 
should be useful to constrain model parameters: 

• The rising and decaying slope of the reverse shock peak. The decaying 
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slope is always in handy since it is the key parameter to identify the 
reverse shock component. It could be used to differentiate the CBM 
profile, e.g. t -2 for ISM and f~ 3 for wind. On the other hand, it 
is also useful to constrain the electron energy distribution index, p r , 
where the subscript r (/) denotes reverse (forward) shock, although 
the constraint is weak, otherwise the decay slope would not be gen¬ 
eral enough for verifying the reverse shock emission. The rising slope 
of the reverse shock is usually missing from the current data, due to 
the limited capability of existing facilities (e.g. slewing speed of the 
dedicated telescopes) and the short-lived nature of the lightcurve rise 
phase. However, once the rising slope becomes available, it is not only 
useful for obtaining the CBM profile, but it is also helpful for test¬ 
ing some proposed extended models, such as the neutron-fed outflow 
model (see details in section [3]) . 

• The reverse shock peaking time is usually related to the shock crossing 
time T x , which is useful to determine the initial physical conditions 
within the GRB ejecta, specifically its Lorentz factor r/ and width A. 
But one needs to keep in mind that the first available observational 
time may not represent the reverse shock peaking time, especially when 
the rising part of the lightcurve is missing. For those cases, only upper 
limits could be made for T x . 

• Based on the standard synchrotron external shock model and assign¬ 
ing reasonable ranges of a set of model parameters, one can constrain 
relevant parameters by fitting the overall observational lightcurve and 
the broad band spectrum, if available. However, in this approach, too 
many unknown free parameters are involved, e.g. the density of CBM, 
the isotropic equivalent kinetic energy of the ejecta, the initial Lorentz 
factor of the ejecta and especially the microphysics parameters in the 
shock region (e e<r ,e ei f, es, r , e e ,/, p r and pf ). Since the observational 
information is usually not adequate to constrain so many parameters, 
some ad hoc assumption are commonly used, for instance, the values 
of the microphysical parameters in the forward and reverse shock re¬ 
gion are assumed the same. It is worth pointing out that the relation 
between e#,r and esj should be treated carefully, since it is useful for 
diagnosing the magnetization degree of the initial outflow. 

• Besides fitting the overall lightcurves, some important parameters such 
as the Lorentz factor and the magnetization degree of the initial out¬ 
flow could also be constrained by working on the “ratios” of the quan- 
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tities for both shocks, especially at T x [581EQ]: 
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where v m , v c and F Ut7n are the typical frequency, cooling frequency and 
the peak flux for synchrotron spectrum, and 
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where 7 C is a critical initial Lorentz factor which divides the thin shell 
and thick shell regimes [58]. This paradigm provides a straightforward 
recipe for directly constraining r] and 1Z B (essentially the magnetiza¬ 
tion degree of the initial outflow) using early optical afterglow data 
only. Moreover, the absolute values of the poorly known model pa¬ 
rameters related to the shock microphysics (e.g. e e , p, etc.) do not 
enter the problem, since they largely cancel out once they are assumed 
to have the same value in both shocks. 


• A morphological analysis of the early optical lightcurves can also pro¬ 
vide direct model constraints. Given a sample of optical lightcurves 
with early detections, one can divide them into different categories 
based on their shapes, then calculate the ratio between each category, 
and find out the right parameter regimes that can reproduce these 
ratios with Monte Carlo simulations [95]. 

• As mentioned above, a time variability within certain modest limits 
in the lightcurve might contain information on some interesting prop¬ 
erties, such as external density fluctuations, inhomogeneous jets, or 
neutron decay signatures, etc. 
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5 Current observational results on reverse shock 


emission 


It has been 15 years since the first prompt optical flash was discovered 
and was interpreted with a reverse shock model (e.g. GRB 990123 [22j f5Tl 
[52]). We have searched the literature since then, finding that 17 GRBs have 
been claimed to have reverse shock signature (3 in the pr e-Swift era). The 
detection rate is much lower than expected. Each of these bursts has been 
interpreted in great detail. In table 1, we collect the burst identifiers and 
their relevant references to the individual studies on those bursts. 


Name 

References 

GRB 990123 

22 52 55 56, 63, 61, 114 115, 116 117, 118 119, 120 121, 122 123, 124 125 

GRB 021004 

25 126 

GRB 021211 

1491 1601 1611 1127111281 

GRB 050525A 

1129111301 

GRB 050904 

I131U132111331 11341 

GRB 06011 IB 

1125111351 

GRB 060117 

ESI] 

GRB 060908 

[S3 

GRB 061126 

1138111391 

GRB 080319B 

11401 11411 [142111431 

GRB 081007 

ESS 

GRB 090102 

1145111461 

GRB 090424 

ESS 

GRB 090902B 

11471 11481 

GRB 091024 

mi 

GRB 110205A 

11501115111152] 

GRB 130427A 

1153111541 11551 11561 


Table 1: GRBs with claimed reverse shock signatures and the corresponding 
references. 

Most recently, a comprehensive statistical analysis of reverse shock emis¬ 
sion in the optical afterglows of GRBs was carried out [157] . Here we briefly 
summarize the results as follows: 

• With stricter criteria, such as requiring redshift measurement, a full 
sample of 10 GRBs with reverse shock signatures was identified: GRBs 
990123, 021004, 021211, 060908, 061126, 080319B, 081007, 090102, 
090424 and 130427A. For five of them, a reverse shock component 
has been firmly confirmed (e.g., GRB 990123 [51], GRB 021211 [25) 
[59], GRB 061126 [1383, GRB 081007 [HI], GRB 130427A [T53llT58] h 
For the remaining five cases, different interpretations (other than the 
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reverse shock emission) can be applicable for the early observational 
results, due to the lack of good early-time photometric coverage. 

• In the sample, GRB 012004 is the only case with a possible Type I 
lightcurve (in which both reverse and forward shock afterglow lightcurve 
peaks are observed) and the other nine cases are all with Type II 
lightcurves (in which the characteristic steep-to-shallow light curve 
evolution is observed). 

• Based on the analytic reverse shock plus forward shock model, the 
physical quantities describing the ejecta and CBM are explored by re¬ 
producing the observed optical lightcurves of the sample with Monte 
Carlo simulations, with the result that the physical properties cover a 
wide parameter space and do not seem to cluster around any prefer¬ 
ential values, which is consistent with previous analyses that concen¬ 
trated on late time forward shock emission 11)1 [20| . 

• It is suggested that GRBs with an identifiable reverse shock compo¬ 
nent show high magnetization parameter i?B = £B,r/^B,f ~ 2 — 10 4 . 
Together with the fact that 9/10 of the cases in the sample belong 
to Type II, the results are in agreement with the mildly magnetized 
baryonic jet model of GRBs m- 


6 Summary and prospects for reverse shock stud¬ 
ies 

Reverse shock emission is a natural prediction of the standard external shock 
GRB afterglow model, and it has been firmly confirmed in a small number 
of cases. Since the reverse shock emission is directly related to the GRB 
outflow itself, investigating the nature of reverse shock emission would lead 
to a better understanding of the intrinsic properties of the GRB ejecta, 
which is essential for constructing a complete picture of the GRB physics. 

A theoretical framework for the behavior of the reverse shock emission 
under various conditions was developed, mostly before the launch of Swift 
(and even before the first relevant discovery of GRB 990123), and expected 
features were discussed for inferring various intrinsic properties of the GRB 
ejecta. Swift was launched, in part, with hopes to make significant progress 
on this specific problem. After a decade of highly successful operation, Swift 
indeed has collected a good sample of early afterglow lightcurves to allow 
detailed studies of GRB reverse shocks. While the size of the sample is 
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still limited, nonetheless it appears that the number of bursts with con¬ 
firmed reverse shock components is much lower than the expectation from 
the standard model. 

The mismatch between this theoretical expectation and the observations 
could be intrinsic or it could be systematically biased due to the limitations 
of current ground-based observational facilities. If it is intrinsic, the origin 
of the suppression of the reverse shock emission for most of GRBs would 
shed new light on the composition problem of GRB jets, e.g., most of the 
jets might be highly magnetized. 

Based on current observational results, more reliable results could also 
be achieved by including more broadband or more specialized information 
instead of just photometric or spectroscopic optical data. For instance, one 
could use early radio data or (sub)mm data [ 159 . (I60j to search for re¬ 
verse shock emission signatures [1531 [161) ; one could identify the reverse 
shock components and diagnose the structure of the magnetic fields in GRB 
ejecta via the detection of early time optical polarization [108 , [109 ); one 
could estimate the magnetization degree of the GRB jets by comprehen¬ 
sive considering the 7 -ray spectrum [162) , the early optical lightcurve type 
and special X-ray afterglow features, such as the X-ray plateau due to late 
magnetic energy injection m- 

At this point, the main problem is that there is still a large fraction 
of GRBs lacking early optical observations, and a more complete sample 
is required for firmer conclusions. Some upcoming facilities may help with 
this issue, such as the Chinese-French mission SVOM [163 ) and especially 
its key element, the Ground Wide Angle Cameras (GWACs). The GWACs 
is an array of wide field of view (about 8000 deg 2 , with a sensitivity of about 
15 magnitudes at 5 s) optical cameras operating in the optical domain. It 
will monitor continuously the field covered by the SVOM 7 -ray detector 
ECLAIRs, in order to observe the visible emission of more than 20% of the 
events, at least 5 minutes before and 15 minutes after the GRB trigger. This 
and other ground-based facilities may key in making further progress in this 
field. 
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